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Introduction 


- In the theory of proton-proton scattering Coulomb wave functions play a 
fundamental réle. For this and similar reasons it seems to be justified to 
examine methods for their numerical computation. Well-known series expan- 
sions have been constructed by Sex [1] and Yost, WHEELER and Breir [2]. 
Further Mathematical Tables Project! (now under the direction of Bureau of 
Standards, Washington, D.C.) is preparing a table of the regular and irregular 
S-wave functions. With this table it will be possible to obtain any reasonable 

- function value with about 6 digits accuracy. 

It is the aim of this paper to give a representation of the Coulomb wave 
functions with real integrals and to show that they can actually be used for 
numerical purposes. The formulas might be useful especially in the case when 
just one or a few function values are needed. 


1. The equation which determines the functions in question can be written 


du (1 2a Uae \u=0. (1) 


de? 0 0 


A fundamental system is defined with the help of the following boundary con- 
ditions: 


; l 
u(0)=0; U(O)o+o0 ~ sin (. Sar (0. «)) (2 a) 


l 
0 (Q)yree ~ C08 (. —in+nle, a). (2 b) 


n(o,a) can be determined from (2a); then v is defined by (2 b) with the 
same 7 (0, a). 


1 IT am indebted to Dr. A. N. Lowan, Math. Tables Project, New York, for placing at 
my disposal the manuscript of the regular Coulomb wave functions table [3]. 
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Putting * 
u= ele 1G (0) (3) 


in (1) we obtain: 


g (G+ 246") + 21+ 1)@ + 1G) —2aG=0. (4) 


On this equation the usual Laplace-transformation is applied: 


G= [4 (B)- e-dB. (5) 


Observing that 9 e?? = 75 e?° we obtain after a partial integration and putting | 


the integrated part and the remaining integrand both equal zero: 
[(B? + 22 B)- > (B)- e7e]2 =0 (6 a) 


#' (6) _ 21(8 + 1) —2a 
(8) B+ 2ap (S2) 


where the limits s; and sg have to be chosen conveniently. Thus 


# (B) = FP + 21 pylertearete 340, (7) 


The integrated part vanishes when 8 = 0, 8 =— 212 and B =— oo (since 9 > 0). 
The first two points are singular, but nevertheless it is possible to choose the 
integration path from one of the points to the other one. (This could be seen 
by taking an integration curve, encircling the two points, and then contracting 
the curve to the double line between — 27 and 0.) 

Putting 6 + 1 =7t and making some simple reductions we deduce: 


1 
1 3 
eet tte)? | (— aF o0e (e¢— a tog alee (8) 


0 


Ue 


f(a) 1s an amplitude factor which will be determined later. When a =O, 
u must go over into 


|/ ture). 
1 


ge | (1 — #7)! - cos gt-de= 2-1! |/ Sy Jt (0) 


Now we have: 


(cf. [4] p. 46), from which we see that /(0)=1. 
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Analogously if we choose the limits — co and 0, we see that we must sub- 
tract 7: wu in order to get the right function v: 


= 7 J a(a)| { (1 ie ry e~ett+2aarctgt Jy — 
; 0 
1 
— [=e sin (o¢— a tog 1) ae}. (9) 
0 
jco pp =—- — eS eS SK Se See ee ee OO —— >} fay ani ie~) 

‘ | 

| | 

| 

| 

\ 

0 Pp 


Vannes, Ue 


2. The asymptotic behaviour for large values of @ can be found by writing: 


Q 
eps ; FR aye ayes (CES) 1 
LS eg a f (a) RP. | (0 Nee ig aa. 
0 


- We replace the integration path along the real axis by 


(0,700) + (00, 9.4 100) + (9 + 100, g) 
(see fig. 1). |; 
Because of the factor e'’* the integral between 1co and @ + 700 drops out 


| and we obtain, putting c=7& and «x =0 + 76 respectively: 


sie AN 


ao 


2 0)—! : —i+2e aretg ~ 
| w= s(OH Rei f et + Hye aide 


— [ 2g gteee Co + ieyte(—igie-iae | 
0 


Re P) = real part. 
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The first integral, being purely imaginary, disappears and we are left with: 


saa ~ : l—ia 2 
ie a RP. ‘ ei (o-a@ log 2 9—34 U7) [ (1 ab | éltta ed :| 5 
; 0 


v 


For large values of @ the factor 1 + may be replaced by unity, and the: 
integral reduces to (ia +1+ 1). Thus: | 


~ 


f(a aetat| al /(l+a 2)(4+ a®)--- (2 + a) na 
/ 


San ‘sin (0 — a log 20 — 41a + m)) 


where 4, = arg (ta + 1+ 1). Now f(a) should be chosen to make the ampli-: 
tude unity, and then we finally get the following formula: 


w= 2- Vs 1 iS git. 
2ma(1 + a?)-:: (2 + a?) 


1 


[ey e0s (gta tog 74) ae. (10) 


In an analogous way we get, finding g(a) =/ (a): 


ioe) 


- il e72me 
— 9-1 Vy os 41 || 2\l p-et+2earetgt, Jy 
Z | 2na(1 + a2) --- (2 + a) “ naa a! 


ia) ae]: (11) 


— | (1 — #7)! - sin (ora ieee 


0 


1 


The formulas (10) and (11) are, of course, essentially identical with the well- 
known integral representations, given e.g. in Morr—Massry [5]. It does not 
seem to have been pointed out before, however, that numerical computation 
actually can be performed on the basis of the formulas above. We will return 
to this later and also give some examples. 


3. For reference we also give the behaviour at the origin. Then we need 
the integral 


, erik 
a | (1 —#)'- cos (« log arma | iy (; a Vat. 
A jicscoy =e 
: i 
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After the substitution ¢ = 2a2—1 the integral transforms to: 


1 
ef ri+1+ia)-Pii+1—ia) 
T= 22! Itia l-ia 20. — 
= les ie i Ie, Tere 3 


0 
Oa ee ae ei. 
(21+ 1)! shza 


The remaining integrals needed are elementary. 
Thus we have the following approximate formulas for small values of 0: 


£222 a(lto)-- (2+) 1, 
/ 3 o) 
ae 90+ 1)! V2 e2te 2= jl Q (12 a) 
£9 e27e |] : 
eee Hp ») 
x ss an Qna(l + a*)--- (P+ a) S te 


(12 a) and (12 b) can easily be obtained from the series expansions (cf. e.g. 
[5], p. 53). Now we also find the equalities: 


ij du 1+] 

Ar ie ae the ol 

- avy l 

lm wg = aT eT a) 
UI VI 1 


4. In the applications the case 1 = 0 will be of special interest. The func- 
tion v, however, contains an integral over an infinite interval, so we make 
the following transformation: = e~°' and replace & by t. Then we get: 


1 


1 — e274 : 1+t 
ag ———"9 | cos @ — a log Te ) dt (14 a) 


0 


5S . log 1/t 
— e274 2a arctg f Oh sees 
Vy) = is s = ay E eaey nae 0 sin (0 — a log +) -dt. (14b) 


270 
0 


Both integrands are’ finite, and even the singularity in ¢= 1 will cause no 
trouble at the numerical computation. The integrand of (14 a) is presented 
m fig. 2. 

A few function values have been computed for a= } and some of them 
could be compared with the corresponding values from the table [3]. The 
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Fig. 2. 


numerical integration was performed by Simpson’s rule and rather large inter- 
vals. There would be but little difficulty to increase the accuracy consider- 
ably. The table below gives the result. 


a=. 
i 
0 uo (comp.) uo (tab.) Diff. % u, (comp.) 
0.1 0.039511 0.039519 0.020 0.00143785 
0.6 0.284800 0.284742 0.020 0.056582 
1.0 0.516583 0.516601 0.003 0.162767 
2.0 1.021184 1.021120 0.006 0.618514 


Institutionen for mekanik och matematisk fysik, Lund. 
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